Search for: All records

Documenting the processes and timescales of magma formation and diversification and defining the locations, shapes, volumes, and phase states of magma storage and transport zones rely on data produced by novel analytical techniques and state-of-the-art experimental methods. Computational modeling effectively links these critical tools. The Magma Chamber Simulator (MCS) is an internally consistent thermodynamic open system model that uses experimental constraints from rhyolite MELTS (Gualda et al. 2012, Ghiorso & Gualda, 2015) to compute paths of open system magmas that evolve via processes including crystallization, magma mixing, cumulate/mush entrainment, and host-rock assimilation. MCS results yield elemental, isotope, mass, and thermal characteristics of melt ± crystals ± volatiles in "resident" magma, crustal wallrock (melt and solids), recharge magma, and entrained material. To model the petrochemical evolution of igneous rocks related by open-system processes, one typically runs 200+ models that vary initial compositions, pressures, and temperatures of magma, wallrock, etc. Comparison of model results with whole rock, mineral, and melt inclusion chemical data and other constraints (e.g., thermobarometry) yield interpretations about igneous processes at a range of scales—from how crust forms and evolves to processes responsible for in situ geochemical records of crystals—and allows assessment of epistemic and aleatoric uncertainties. Two examples of computational studies will illustrate MCS's utility and flexibility. (1) Modeling of historical basalts at Mt. Etna (Italy) provides evidence for variable degrees of melting of metasomatized mantle, followed by magma recharge and assimilation of partial melts of carbonate-flysch crust. (2) MCS models reproduce whole rock and mineral data of plagioclase-rich basalts at Steens Mountain (USA) through entrainment of gabbroic mush that likely formed in early stages of Columbia River Basalt magmatism. To enhance understanding of trans-lithospheric magma systems, future work on MCS will prioritize (i) building a post-processing environment that utilizes select statistical methods to inform "best-fit" models and to quantitatively assess uncertainty, and (ii) increasing modeling efficiency by adding automated modeling capabilities. 

Several times in the past 60 ka, Mount Etna has erupted lavas with variable alkaline character. The most recent chemical excursion began in 1971, accompanied by an increase in explosivity and eruption frequency. The origin of the alkaline signature remains enigmatic, with endmember hypotheses involving dominant contributions from mantle vs. crust. For lavas that erupted between 1329 and 2016, we used thermodynamic modeling to test if post-1971 anomalous alkalinity is dominated by mantle processes. First, we assessed mantle melting conditions required to reproduce the chemistry of potential parental magmas. Second, we examined the differentiation of partial melts as they underwent closed-system crustal storage and ascent. The mantle melting conditions explored via ~300 models include source compositions (peridotite + pyroxenite ± metasomatic phases), pressure, extent of melting, and fO2. Best-fit models that reproduce the major element chemistry of volcanics interpreted as parental magmas involve 20 to 30% melting of a peridotite-pyroxenite mantle source that contained phlogopite, pargasite, and CO2 (~1 wt.%) between ~1-1.5 GPa (~30-45 km) along the QFM+1 buffer. Subsequent isentropic decompression (adiabatic and reversible) of mantle partial melts + crystallization at shallower pressures (0.8-0.2 GPa) were modeled to test the effects of closed system ascent and storage. Isentropic decompression models yield no crystallization although temperature decreases ~3°C/0.1 GPa. Decompression + closed system crystallization fail to replicate post-1971 glass samples and do not explain observed post-1971 alkali enrichments. We conclude that partial melting of a metasomatized source produced Etna parental magmas, but closed system crustal ascent and storage cannot fully account for alkali enrichment highlighted in post-1971 products at Etna. Open system modeling suggests that assimilation and crystallization (e.g., Takach et al. 2024) play a critical role, and ongoing modeling is testing the contributions of recharge (magma replenishment) and entrainment of previously formed mushes to the high alkalinity excursions.